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Abstract: A new closed-loop drive scheme which decouples the phase and the gain of the 
closed-loop driving system was designed in a Silicon Micro-Gyroscope (SMG). We 
deduce the system model of closed-loop driving and use stochastic averaging to obtain an 
approximate "slow" system that clarifies the effect of thermal noise. The effects of 
mechanical-thermal noise on the driving performance of the SMG, including the noise 
spectral density of the driving amplitude and frequency, are derived. By calculating and 
comparing the noise amplitude due to thermal noise both in the opened-loop driving and in 
the closed-loop driving, we find that the closed-loop driving does not reduce the RMS 
noise amplitude. We observe that the RMS noise frequency can be reduced by increasing 
the quality factor and the drive amplitude in the closed-loop driving system. The 
experiment and simulation validate the feasibility of closed-loop driving and confirm the 
validity of the averaged equation and its stablility criterion. The experiment and simulation 
results indicate the electrical noise of closed-loop driving circuitry is bigger than the 
mechanical-thermal noise and as the driving mass decreases, the mechanical-thermal noise 
may get bigger than the electrical noise of the closed-loop driving circuitry. 

Keywords: Silicon Micro-Gyroscope (SMG); drive-mode; thermal noise; stochastic 
averaging 
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1. Introduction 

The Silicon Micro-Gyroscope (SMG) is an important MEMS inertia sensor with a broad application 
in the national economy and defense fields [1-7]. These SMGs usually measure the proof-mass 
displacement by capacitive methods, but under normal atmospheric pressure, the minute moving 
proof-masses are especially susceptible to mechanical noise resulting from molecular agitation. 
Although accuracy is usually limited by electrical noise and systematic errors, mechanical thermal 
noise provides a theoretical lower limit for random errors [8-11]. Thus, a proper accounting of thermal 
noise is essential for the development of higher accuracy tactical and inertial grade gyroscopes. 

The effects of mechanical-thermal noise on the sense-mode have been presented in the literature [8- 
11], but discussions of the effects of mechanical-thermal noise on drive-mode can hardly be found in 
the current literature. In this paper the effects of mechanical thermal noise on the driving performance 
of the SMG are mainly derived. Only the influence of the mechanical thermal noise is considered, 
while the electrical noise, sampling and quantization error, and distortion due to filtering are not 
considered. Meanwhile, we assume all the other processes run in an ideal manner. In this paper, a 
stochastic averaging approach is used to take account of the effects of closed-loop drive. The effect of 
mechanical thermal noise on drive-mode is discussed, then stochastic averaging is used to develop a 
model for the "slow" dynamics which represent the driving amplitude and frequency of the SMG. Both 
the steady-state and transient response of the model are obtained by stochastic averaging. The spectral 
density of the random error due to thermal noise on drive-mode is also derived. 

2. Working Principle 

As shown schematically in Figure 1(A), the micro-gyroscope consists of two silicon frames (outer- 
frame and inner- frame); the outer- frame is anchored on a glass substrate by six outer support beams 
and is connected with the inner- frame through four inner support beams. The outer-frame and the fixed 
interdigitated drive electrode on the glass substrate form the drive capacitors. The alternating drive 
force of the out- frame along the x-axis is generated through applying alternating current (AC) voltage 
with direct current (DC) bias voltage to the fixed drive electrode. Since the stiffness of the inner 
support beam along the x-axis (K x i » K x ) is very large, the outer- frame and the inner- frame are driven 
together to vibrate along the x-axis by the alternating drive force, which causes the alternating 
capacitance between the outer-frame and fixed drive-sense electrode. We can capture the drive 
displacement by detecting the alternating capacitance. When the rotation rate along the z-axis is input, 
according to the Coriolis effect, the Coriolis force along the y-axis will be loaded on both the outer- 
frame and the inner- frame. Because the stiffness of the outer support beam along the y-axis (K yo » K y ) 
is very large, only the inner-frame is driven to vibrate along the y-axis by the Coriolis force, which 
induces the alternating capacitance between the inner-frame and fixed sense electrode. We can obtain 
the rotation rate along the z-axis by detecting the alternating capacitance. 



Sensors 2009, 9 



3359 



Figure 1. (A) The frame of the SMG. (B) The simple model of SMG. (C) The picture of 
the processed SMG. 

(A) (B) (C) 




The simplified motion equations of SMG are described by: 

m x x + R x x + K x x = F e +n (1) 
m y y + R y y + K y y = -2m y Qx (2) 

where x and y are separately the drive axis displacements and sense axis displacements in meters, £1 
the rotation rate along the z-axis in radians/second, m x (m x =mi+m2) and m y (m y =m2) the drive proof 
mass and the sense proof mass in kilograms, R x and R y the damping in Newtons/meter/second, K x and 
K y the stiffness in Newtons/meter, and-2m x Qx the denote of the Coriolis force. F e (F e =FdsincOdt) is the 

electrostatic force used to maintain the drive-mode vibration at a specified amplitude in terms of 
displacement, and at a resonant frequency of the drive-mode. Mechanical thermal noise on the drive 
axis is represented by the random force n(t), in units of force. 

Ignoring the influence of the random force n(t), the drive axis displacements and sense axis 
displacements in the steady state are described by: 

x « A x sm(co d t + (p) (3) 
y*A y sm(co d t + (p + y) (4) 

where A x = , , ; (P = ~tg ( — ; 2 2 ) ; A y * — * 



If 2 2 \2 C0 2 nx CO d 2 Qx^nx-Vd) / 2 2 X2 (0*0)/ 



y = ^-i ( ^KlA1) ; co nx =(K x /m x ) (1/2) ; co ny =(K y /m y ) (1/2) ; Q x =m x co nx /R x , Q y =m y co ny /R y 



When C0d=c0nx = 03 n y, the maximum drive axis displacements and sense axis displacements are 
described by: 

4) = ^sin( ffl /~) (5) 

-2F.Q Q Q n 
y(0= f My ^{co d t~) (6) 
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3. Mechanical Thermal Noise 

Consider the damped harmonic oscillator: 

m x x + R x x + K x x = n (7) 

The presence of damping in the system suggests that any oscillation would continue to decrease in 
amplitude forever. Inclusion of the fluctuating force n(t) prevents the system temperature from 
dropping below that of the system's surroundings. The damper provides a path for energy to leave the 
mass-spring system. This is the essence of the Fluctuation-Dissipation Theorem. According to 
Equipartition, if any collection of energy storage mode is in thermal equilibrium, then each mode will 
have an average energy equal to (l/2)k B T where k B is Boltzmann's constant(1.38><10" J/K) and T is 
the absolute temperature in degrees Kelvin. A mode of energy storage is one in which the energy is 
proportional to the square of some coordinate ; e.g., kinetic and spring potential. 

When this system is in thermal equilibrium, the probability distribution of x and x is given by 
equation (8): [12] 

E(x,x) 

p(x, x) = const ant x e keT (8) 
For the oscillator, the energy is the sum of the kinetic and spring potential energy: 

E(x,x) = ^m x x 2 +^K x x 2 (9) 

From here, the equipartition theorem can be derived, namely that the mean energy in any energy 
storage mode is equal to(l/2)k B T. Thus: 

< - m Y x >=< - Kx 2 >= — k R T (1 0) 

2 2 2 

where < >denotes an ensemble average. The form of the distribution p(x, k) indicates that x and x are 
independent, Gaussian, and have zero mean. Since this holds for all values of m x and K x , the thermal 
noise n(t) must be a white Gaussian noise with two sided spectral density [10]: 

S n (co) = 2k B TR x N 2 /Hz -oo<«<+oo (11) 



The spectral density of drive displacements due to thermal noise is: 

\_ 

{col- co 2 ) 



SJ co) =— * — S„( co) ( 1 2} 

2 „2\ 2 + CO nx CO 



rn 2 



Q 



J 



So the noise power spectrum is: 

<x 2 >=—\ S x ico)dco = ^- (13) 



2k i - K 
The RMS noise displacements due to thermal noise is: 
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According to Equation (5) and Equation (14) the signal-to-noise ratio of drive-mode is: 

(S/N) x =x/x n =^&= (15) 

The above equation indicates that we can improve the signal-to-noise ratio by increasing the 
quality Q x and driving force amplitude Fd, or by reducing the stiffness and temperature. 

4. Closed-loop Driving and System Model 

As is known in the art of Coriolis force sensors, in order to achieve an acceptable response from the 
sensor, the proof mass vibration of the drive-mode should have a frequency at, or close to, the resonant 
frequency of the proof mass. At the same time, in order to improve the entire performance of the SMG, 
a high stability of the driving frequency and the amplitude of the drive-mode are needed. To satisfy 
those demands, the closed-loop driving of the drive-mode must be achieved. To this end, the drive 
signal has a frequency equal to the resonant frequency of the proof mass. However, parasitic 
capacitances between the drive electrode and the drive-sense electrode can cause significant errors. 
That is, when the drive signal capacitively couples into the drive-sense electrode, the accuracy of 
amplitude control by the feedback circuit is degraded and the harmonic frequency of the closed-loop 
system departs from the resonant frequency of the proof mass, resulting in less than optimum sensor 
performance, so we must eliminate the capacitive coupling. Various techniques are generally utilized 
in an effort to reduce capacitive coupling. In this paper, such a technique is utilized as follows: the 
drive electrode is arranged on the left, the drive-sense electrode is arranged on the right and the anchor 
of the SMG is connected with the ground or the virtual ground, which is shown in Figure 1(A). By 
separating the drive-sense interface (drive-sense electrode) from the interference source (drive voltage), 
we can reduce these capacitive couplings. Besides, the SMG studied in this paper is executed in 
vacuum encapsulation, with a working pressure under 10" 1 Torr and the quality factor of the drive- 
mode above 2,500, which can also reduce these capacitive coupling by decreasing the drive voltage. 
The modulation-demodulation method through applying high-frequency carrier to the proof mass can 
also reduce these capacitive couplings. 

First, we need to extract the resonance signal of the drive-mode. The simplified interface circuitry is 
shown in Figure 1 . Figure 2 show the details of the interface circuit and the equivalent circuit. Here in 
Figure 2, C(t) is the alternating drive-sense capacitance, C(t)=Co+AC, Co is constant capacitance. The 
part of signal sense can be equivalent to a current supply I(t) and a internal resistance Co [see Figure 
2(B)], where: 

m _dQ_ d(y,C(t)) 

In Equation (16), the differential of driving capacitance to displacement K xc (K xc =— ) is a 

dx 

constant relating to the structure. The capacitance Co is very minute and generally has hundreds of fF, 
thus the impedance of Co is very large and we can ignore the influence of the impedance Co. The 
resistance Ri generally has a few Mf2, the capacitance Ci hundreds of nF, co/ 2JI a few kHz (co~co nx ), 
so Ri»l/coCi, the output voltage P(t) is described by: 
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P(t) = I 1 (t)R = I(t)R 



R 



R { +1/ jcoC x 



* V.RK„ x 



(17) 



Figure 2. (A) The interface circuit of drive-sense signal. (B) The equivalent circuit. 
(A) (B) 
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Figure 3 is the frame of the closed-loop driving. In this figure s is a complex variable. V s , V re f and 
V sup are the direct current biases. L and J are the zero and the pole of the integrator separately. G is the 
gain of the integrator. In order to improve the precision and the stability of the closed-loop driving, the 
Q-factor of the drive-mode should be increased (the SMG is executed in vacuum encapsulation), while 
the well closed-loop control should be achieved. The closed-loop control must meet such two 
conditions: 1. The phase of the whole loop Q=2rm (n is an integer); 2. The gain of the whole loop A>1. 



Figure 3. The frame of the closed-loop driving. 
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A new closed-loop drive scheme which decouples the phase and the gain of the closed-loop driving 
system is adopted in the SMG, so that the phase and the gain can be optimized, respectively. The gain 
of the whole closed-loop system is controlled by the branch circuit above, and the phase is controlled 
by the branch circuit below. These two branch circuits respectively fulfill the two conditions of closed- 
loop control, adjusting and optimizing the closed-loop parameter separately. The "voltage comparator" 
is the key component of the closed-loop driving. The output of the "voltage comparator", with an 
invariable output amplitude, only reserves the phase information of the input signal, so the phase 
conditions of the closed-loop are isolated from the gain conditions. Except the drive mode of SMG and 
"voltage comparator", with the suppose that the phase of the other parts is fixed, when the phase of the 
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"voltage comparator" is changed, the vibrating frequency of the closed-loop system will depart from 
the resonant frequency of the proof mass, so the phase condition controls the frequency of the driving 
displacement. It is obvious that the gain branch controls the amplitude of the driving displacement. 
According to Equations (34) and (35) (See Section 5), ignoring the influence of the random processes 
ni(t) and n 2 (t) ,we can know that average amplitude a and average phased are decoupled, so the phase 

branch and the gain branch can be optimized respectively and the new closed-loop drive scheme is 
succeeded. 

To make sure that the harmonic frequency of the closed-loop system equals, or gets close to the 
resonant frequency of the proof mass, that is to say, the phase of the whole loop 0 = 2n7i(n is an 
integer). When co d =co n x, the phase-shift of the drive displacement x(t) comparing to the drive force F e 
is -7i/2(See Equation (5)), the phase-shift of the output voltage of pre-amplifier P(t) comparing to drive 
displacement x(t) is -nil and the phase-shift of the output voltage of the "voltage comparator" V cp (t) 
comparing to output voltage of preamplifier P(t) is -7i. The other parts, which in fact all have tiny 
phase errors, have no phase-shifts, so the closed-loop control meets the phase conditions. In this way, 
it is hoped that the above closed-loop phase errors should be as tiny as possible. Various techniques 
are generally utilized in an effort to reduce closed-loop phase error, or drift, in servo circuits, such as 
amplifier circuits utilizing an operational amplifier. One such technique includes the addition of one or 
more zeros (i.e., a lead filter) in cascade with the open-loop gain of the operational amplifier in order 
to flatten the open-loop gain over a portion of the frequency band, generally resulting in only moderate 
closed-loop error reduction and also compromising stability. Another technique for reducing gain and 
phase errors is to increase the gain-bandwidth product associated with the operational amplifier. 
However, use of this technique is limited by the gain-bandwidth product of commercially available 
operational amplifiers as well as by the acceptable increased power dissipation associated with higher 
performance operational amplifiers. However, a Phase-Corrected Amplifier Circuit can be used to 
remove the closed- loop phase error [13]. An amplifier circuit having a bandpass circuit in cascade with 
the forward loop gain is provided, with the bandpass circuit having a transfer function approximating 
one plus a bandpass characteristic, the passband of which corresponds to the information band. This 
arrangement increases the open-loop gain of the amplifier circuit around the information frequency 
without affecting the open-loop gain at DC and crossover so as to reduce phase and gain errors around 
the information frequency. 

In Figure 1, the drive capacitance and the electrical potential energy stored in the capacitance is 
described as: 



where So is the permittivity, h the thickness of the comb fingers, xo the overlap length of the fingers, 
and d the width of the gap between fingers. According to Figure 3, the electrostatic drive force is 
described by: 



Ql(0 = *0 



h(x 0 + x{t)) 
d 



(18) 



U=±C dl (t)V 2 



(19) 



1 dC dl (t) 



1 dC dl (t) 



[(V sup -z(t))V cp (t) + (V sup + z(t))] 



2 dx 



2 Dx 
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= r^r^- +z{t)f + ih (v * lp ~ z ^ v ^ + \^ 1( y^ - ^))%(o a (20) 

3C 

where — —is the differential of driving capacitance to displacement, the z(t) is a DC and V cp (t) is a 

dx 

AC with the frequency near the resonant frequency co nx . The first term of Equation (20) is a constant 
and does not contribute to the oscillating driving force. Since the amplitude of the AC driving voltage 
is chosen to be much smaller than the bias voltage, the third term of Equation (20) is much smaller 
than the second term and can be neglected. Therefore, under proper bias, the driving force is 
approximately proportional to the AC driving voltage. 

According to Figure 3, the drive-mode of the SMG is modeled as a second-order spring-mass- 
damper system with a dynamic behavior described by: 

x + ^x + ^x = F>) + — (21) 

DC 

where F n (u) is electrostatic drive force (F (u) = F e I m x « — — 0^L ~z(t) 2 )V (t)/m ), and suppose 

q x 

z(t)>0. 

Simplifying the integrator (See Figure 3) into the basic integral function, so: 



-G\[(V ref -\P(t)\)dt = z(t) 



i.e.: 

z(t) = -G(V ref -V s R\K xc \\x\) (22) 

where | x | is modulus of the driving velocity and: 

P(t) = V s RK xc x 

So the output of the "voltage comparator" is: 

V (0 = --^- = - KJt (23) 
cA) \P(t)\ \K xc \\x\ 1 ^ 

In summary, the entire closed-loop driving system, shown in Figure 3, are described by Equations 
(24)~(25): 

x + R x + at x x = — (24) 
m x 

z(t) = -G(V ref -V s R\K xc \\x\) (25) 



where Re q (R eq = + (V sup - z(t) ) — — ) is equivalent damper. Equations (24) and (25) 

Q x dx I K xc \\x\m x 

show the closed-loop system model. With the ignorance of the influence of the random process n(t), 
Equation (24) shows a free spring-mass-damper oscillating system and Equation (25) shows the 
control principle of the integrator. The closed-loop system reduces the system damper Re q through 
adjusting the output of the integrator z(t). When the equivalent damper R eq is bigger than the 
zero(e.g.Req>0), the gain control loop reduces the output of the integrator z(t) and then the equivalent 
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DC 

damper R eq will decrease near the zero (R eq ~0, — — is inverse-phase with K xc ). When the equivalent 

dx 

damper R eq is smaller than the zero (e.g.R eq <0), the gain control loop enhances the output of the 
integrator z(t) and thus, the equivalent damper will increase near the zero (Re q ~0). So the closed-loop 
system shows approximately an undamped-free vibration with the invariable amplitude and 
frequency(e.g. the resonant frequency co nx ). However, the random process n(t) will influence the 
stability of the driving amplitude and frequency. 

5. Stochastic Averaging of the Gyroscope Dynamics 

After the vacuum encapsulation, the quality factor of the drive-mode Q x becomes very big. The 
drive-mode of the SMG can be equivalent to a band-pass filter. Only the displacement with the 
resonant frequency co nx can be magnified and the other harmonics are attenuated greatly. So the 
transient displacement can be simplified into a pure sin component. In order to analyze the transient 
behavior of the system, the driving displacement of the SMG is defined as [14]: 

x = a(t) cos(co nx t + <f>(t)) (26) 

where a(t) is the amplitude, and </>(t) is the phase of the driving displacement signal respectively. 
Differentiating Equation (26) with respect to time gives the velocity: 

x = -aco nx sm(a> nx t + (/)) + a cos(co nx t + (f))-a(j) sm(a> nx t + (/)) 

One of the equations used to determine a and (f) is obtained by assuming that the sum of the last 
two terms can be set equal to zero: 

a cos(a> nx t + (fi)-a^ sm(co nx t + 0) = 0 (27) 

Thus, the velocity equation becomes: 

x = -aa> ru sm(a) rtt t + 0) (28) 

The acceleration is obtained by differentiating Equation (28) with respect to time: 

x = -aco nx sin(^/ + (/))- aco nx {co nx + j>) cos(co nx t + (/)) (29) 

Substituting Equations (26), (28) and (29) into Equation (24) yields: 

co 2 • K asiwico t + d>) n /^rw 

_ i a O}nx+a ^)sm(a } J + 0)-aa } Jcos(a } J + 0) = ^(V^ + zmV sw -z(t)) f V > + — (30) 

Q ox K \\asm(a) t + 6)\m m. 

I xc 1 1 v nx r / \ x X 

Substituting Equation (28) into Equation (25), and suppose a ^0, according to Equation (27) and 
Equation (30), separate a , (j> yields: 



<t> = ~ 7f sm(»„; + (/>) cos(coj + (/>)- — - — — — p p - 

Q x dx I K xc || a sm(a)J + </>) \ m x a) m 
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-cos(G)J + 0) (32) 



am co 



z{t) = -G(V ref -V S R | K xc || a«„ x sin(^ + 0 |) (33) 

It should be noted that Equations (31)~(33) are the exact differential equations describing the 
evolution of the amplitude and phase of the driving displacement, as well as that of the output states of 
the integrator. However, these equations are difficult to analyze because they are substantially non- 
autonomous. It is evident that instantaneous phase co nx t evolves much faster than the other variables, 
such as a, (f) , z and functions sin(co nx t+^) and cos(co nx t+^) are almost periodic. Within a period of 

these sinusoidal functions, variables other than co nx t change very little. Hence, it is possible to apply 
the averaging method to the non-autonomous system described by Equations (3 1)— (33) and 
approximate it by an autonomous system [14-15]. 

As pointed out above, instantaneous phase co nx t is regarded as an independent variable and the 
differential equations Equations (31)~(33) are averaged, with respect to co nx t, over the interval [-71, n]. 
The averaged autonomous equations are: 

t = a co nx 8C dl 2K xc (V suv + I)(V sup - z ) | n, ^ 
2Q x dx n\K xc \m x co nx Sm x co w 



W 

yJ2am x co nx 

U^- W-W ) (36) 

71 

where the bars denote averaged variables. The random processes ni(t) and n2(t) are independent white 
noise with the same intensity as n (t) [10]. Equations (34)~(36) describe approximately how the 
displacement amplitude and phase evolve with time. According to Equations (34) and (35), it is 
evident that average amplitude a and average phase </> are decoupled - when one of them is changed, 

the other does not change, so the phase branch and the gain branch can be optimized respectively. 

Ignoring the influence of the random processes ni(t) and n 2 (t), the equilibrium of the averaged 
system described by Equations (34)~(36) is: 

^= _;;;r, (37) 

ec (38) 

ox 

When the power is switched on, the system stabilizes finally in equilibrium. According to Equation 
(37), it is evident that the equilibrium of the average amplitude a o is independent from the quality 

factor Q x , so the change of quality factor resulting from the variety of temperature and pressure does 
not impact on the equilibrium of the average amplitude a o , which is very important in the practical 

application. The speciality of this system is obviously different from the literature [16]. 

Ignoring the influence of the random processes ni(t) and n 2 (t), the Jacobian matrix of the nonlinear 
dynamic system of Equations (34)~(36) at equilibrium is: 
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df 



4K., 



xc I x nx 



2Q X 7T I K_ | m„6>_ 3x 

2GVR I A" I <y. 



v 2 + 

sup 



*^X J XC 

ox 



xc i nx 



n 



(39) 



and its characteristic equation is: 



5C rfl 8GV„RK v 



2Q x dx n m 



s xc 



i 



V 2 + 

sup 



SQ x V s R Kxc ^ 
ox 



= 0 



(40) 



where y is the variable of the characteristic equation. All the eigenvalues of the linearized averaged 

system are asymptotically stable, if and only if all coefficients is positive. It is evident Equation (40) 

dC 

satisfies this condition because — — is inverse-phase with K xc . In order to make sure that the square 

dx 

root is in existence, the inequation hereinafter must be satisfied: 

dC 



dl 



dx 



)>0 



i.e.: 



„ ^ -Wsl P QxVsRK xc 8C dl 

re f ~ 2 _ _ ^ refo 



dx 



(41) 



where the V re f is the referring DC voltage used as a reference to the amplitude of the pre-amplifier 
output voltage. The V re f 0 is the criterion voltage. When V re f < V re f 0 , the gain control branch works 
normally. When V re f >V re f 0 , the gain control branch loses the control ability. Because the Equations 
(34)~(36) are extremely complex, the averaged system is linearized in the equilibrium, i.e.: 



a = - 



a„ +- 



2Q X U yfll 



- + - 



sc dl 



m co 

x nx 



x\K xc \m x Q) nx dx 



y2 , ^y r ef m x a nx - 
sup Q C Z.. 



ox 



V2< 



am x^nx 



f _ 2GV S R I K xc | e>„ 

71 

where a n =a-a o , z n = z-z 0 .a n can be expressed in Laplace transform domain as: 

- yf2m x co nx 
A n (s) = - 



2Qx 



(42) 

(43) 
(44) 



(45) 
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SGVRK„^ I "I 



where ^ lyi , f ^/"^n* ; ^(s) is Laplace transform of a n , and at (s) is 

J SUP SQ X V S RK rC "" 



3x 

Laplace transform of the ^ , so the steady- state spectral density for the noise component of the a n is: 

V »> = ^ V (46) 

V J 

So the noise power spectrum is: 



< a] >= — \ S-( co)dco = 



In J 

-00 x 

The RMS noise amplitude due to thermal noise is: 



Comparing Equation (14) with Equation (47), we know the RMS noise amplitude due to thermal 

noise in the opened-loop driving is equal to that in the closed-loop driving, so the closed-loop driving 

does not reduce the RMS noise amplitude due to thermal noise. By the way, the RMS noise amplitude 

is independent from the quality factor Q x . 

According to a = a o + a n ~ a o , Equation (43) can be rewritten as: 

a n =$ = n =^ (48) 



o x nx 



S^co) = _ 2 2 2 V<p) (49) 



where co n is the noise frequency, so the power spectral density of the noise frequency co n is 

I 

Suppose the work bandwidth is feHz, the noise power is: 



<co 2 n >=— f S^co)dcD= 2f ^ BTCOnx (50) 

171 '-lis B n a oKQ x 

The RMS noise frequency due to thermal noise is: 

r^T~ V2 \f B k B Tc °nx ^2 / f B k B T 
® n =4<CD 2 n >=— J ,* nx =— Jb b — (51) 

It is useful to reduce RMS noise frequency by increasing quality factor Q x and drive amplitude a o . 



6. Experiment and Simulation 



6.1. Simulation 

In order to validate the feasibility of closed-loop driving and confirm the validity of the averaged 
equation and its stable criterion, the closed-loop driving in Figure 3 is simulated according to 
parameters shown in Table 1 . Figure 4 is the Matlab simulation frame of the closed-loop driving. 
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Figure 4. The Matlab simulation frame of the closed-loop driving. 
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Table 1. The value of simulation parameter. 



Parameter 


value(unit) 


Qx 


2,500 


G 


999 


K X c 


2.92xl0 -8 (F/m) 


V s 


10 (V) 


R 


10 (MQ) 


Vsup 


5(V) 


Vrefo 


0.596 (V) 


L 


-50 


J 


-55 




1.46xlO" 6 (F/m) 




25120 (rad/s) 


m x 


2.89xl0" 7 (kg) 



Figures 5 and 6 are the simulation curves for closed-loop driving when V ref = 0.5 < V refo 
andV ref = 0.7 > V refo [see Equation (41)]. It can be concluded from Figures 5(C) and 6(C) that the 

system resonant frequency approximately equals to the natural frequency of drive-mode after the 
system is stabilized, which indicates that the whole loop's phase 0=2nji(n is an integer)and the phase 
branch circuit below is controlled accurately. When V ref < V refo , according to Figure 5 and Equation 

(37), the driving displacement goes along with the increase of V ref and the gain branch circuit above is 

controlled accurately. When V ref = V refo , the driving displacement reaches the maximum. When V ref 

increases persistently, that is to say, V ref > V refo , according to Figure 6, the output of the integrator z~0 

and the driving displacement does not go along with the increase of V ref and the gain branch circuit 

above does not work normally. The curve of the transient course of drive displacement along with the 
change of time goes smooth and the overshoot disappeares. Comparing Figure 5(A) to Figure 6(A), 
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when V ref < V refo , we can see that the bigger the V ref is, the shorter the transient time is . When 
V ref >V refo , the above gain branch circuit does not work normally and the closed- loop system loses the 
adjustable ability, with the overshoot fading away. 



Figure 5. The simulation curves of closed-loop driving when y = o.5 < V refo • (A) The curve 

of the drive displacement with time. (B) The curve of the output of integrator z(t) with time. 
(C) The curve of the vibrating frequency of the drive-mode with time. 
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Figure 6. The simulation curves of closed- loop driving when y = 0.7 > y refo . (A) The curve 

of the drive displacement with time. (B) The curve of the output of integrator z(t) with time. 
(C) The curve of the vibrating frequency of the drive-mode with time. 
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Figures 7 and 8 are the simulated closed-loop response and averaged equation simulation curves 
when V ref = 0.3 < V refo and V ref = 0.7 > V refo . The simulation result indicates the envelope curve of the 

closed-loop response is basically the same as the curve of the averaged equation. The tiny differences 
come from the simplification of the integrator and truncation error in simulation. From Figure 8, when 
V ref = 0.7 >V refo , we can see that the envelope curve of the closed-loop response is also basically the 

same as the curve of averaged equation. Meanwhile, z ~0, revealing that the above gain branch circuit 
does not work normally, which is in agreement with the practical situation. Thus, the averaged 
equation can be applied to the situations when V ref > V refo through the limit of the working scale of z(t). 
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Figure 7. Closed-loop response (gray) and averaged equation simulation (black) when 
V ref = 0.3 < V refo . (A) The curve of the drive displacement with time. (B) The curve of the 

output of integrator z(t) with time. 




time(s) 



time(s) 



Figure 8. Closed-loop response (gray) and averaged equation simulation (black) when 
y =oj>v refo . (A) The curve of the drive displacement with time. (B) The curve of the 

output of integrator z(t) with time. 
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In Figure 9, the simulation curves of closed-loop driving between Q x =2,500 and Q x =5,000 are 
compared. From Figure 9(A), we can see that although the quality factor increases, the amplitude of 
the driving displacement, which is in agreement with the conclusion derived from Equation (37), 
doesn't change. That is to say, the displacement amplitude has nothing to do with the quality factor. 
From Figure 9(B), we can see that the augment of quality factor results in the growth of output of the 
integrator z(t), which is accordant with the Equation (38). According to Equation (47) and Equation 
(51), when co nx = 25,120 (rad/s), T=300 K, f B =100 Hz, Q x =2,500, a o = 5 |iim, the influence of thermal 

noise on driving performance in different drive proof masses is shown in Figure 10. From Figure 10, 
we can see that the RMS noise displacement and the RMS noise frequency decrease with the increase 
of drive proof mass. 



Sensors 2009, 9 



3372 



Figure 9. The simulation curve of closed-loop driving in Q x =2,500 (black) and Q x =5,000 
(gray) when y = o.3 < V refo • (A) The curve of the drive displacement with time. (B) The 

curve of the output of integrator z(t) with time. 
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time(s) time(s) 

Figure 10. The influence of thermal noise to driving performance when co nx = 25,120 
(rad/s), T=300 K, f B =100 Hz, Q x =2,500, a 0 = 5 jun. (A) The RMS noise displacements 

with driving mass . (B) The RMS noise frequency with driving mass. 
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6.2. Experiment 

The whole experiment circuit is constructed on the idea of Figure 3 and shown in Figure 1 1 . The 
experiment results show that the SMG achieves closed-loop driving and the drive frequency works 
around the natural frequency of SMG all the time. The vibrating waveshape of the closed loop driving 
signal is shown in Figure 1 1(C). It is practical to find a better working point for SMG by adjusting z(t) 
in Figure 3. The experiment results of amplitude control with temperature are shown in Figure 12. 
According to the Equation (5), when cOd=co n x, the amplitude of driving displacement is proportional 
with the quality factor Q x , so the change of quality factor, due to the variety of temperature and 
pressure, impacts directly on the amplitude of driving displacement in the open-loop driving. But the 
new closed loop driving is immunized to the change of the quality factor [see Equation (37) and 
Equation (47)]. With a temperature increase from -40°C to +40°C, the quality factor of the drive-mode 
Q x decreases about 2.6 times [see Figure 12(A)], while the drive amplitude signal increases by only 
0.083% in the closed loop driving [see Figure 12(B)]. 
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Figure 12. The experiment results of amplitude control with temperature. (A) The curve of 
the quality factor of the drive-mode Q x with temperature. (B) The signal amplitude of 
closed loop driving with temperature. 
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The experiment results of closed-loop driving are shown in Figure 13. Figure 13(A) and Figure 13(B) 

are the frequency drift and amplitude drift of driving signal for 1 h respectively. The experiment 

results show that the standard deviation of drive frequency is 0.0205 Hz, with relative drift 5.0 ppm, 

the standard deviation of the amplitude 0.0165 mV and relative drift 14.7 ppm, respectively. Figure 

13(C) is the frequency spectrum of driving signal. The driving signal is about 100 dB bigger than the 

ground noise, which is accordance with the result of the Figure 13(B). As a result, both stabilities of 

driving frequency and amplitude are rather high and the closed loop control is very successful. 
According to Figure 10, when m x =2.89><10" 7 (kg), a o = 5um, co nx = 25,120 (rad/s), we can find that the 

relative noise of drive frequency is approximately 0.0119 ppm and the relative noise of the amplitude 
is approximately 0.8 ppm, respectively. The noise comparison is shown in Table 2. So, in this 
situation, the electrical noise of closed-loop driving circuitry is bigger than the mechanical-thermal 
noise. With the driving mass decreasing, the mechanical-thermal noise may get bigger than the 
electrical noise of closed- loop driving circuitry (see Figure 10). 
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Figure 13. The experiment results of closed-loop driving .(A) Frequency drift of driving 
signal for 1 h. (B) Amplitude drift of driving signal for lh. (C) Frequency spectrum of 
driving signal. 

(C) 
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Table 2. The noise comparison. 



Simulation thermal noise Measure noise(Driving signal) 

m x =2.89xl0" 7 (kg), a o = 5pm, co nx = 25,120 (rad/s) 

Absolute value Relative value Absolute value Relative value 

Thermal noise 

0.0000476 Hz 0.0119ppm Noise frequency 0.0205 Hz 5.0 ppm 
frequency FF H y FF 

Thermal noise 1? 

t 4xl0" 12 m 0.8 ppm Noise amplitude 0.0165 mV 14.7 ppm 
amplitude 

7. Discussion and Conclusions 



The mechanical thermal noise on drive-mode is discussed, and then stochastic averaging is used to 
develop a model for the "slow" dynamics that represent the driving amplitude and frequency of the 
SMG. Both the steady-state and transient response of the model are obtained by stochastic averaging. 
The spectral density of the random error due to thermal noise on drive-mode is also derived. By 
calculating and comparing the RMS noise amplitude due to thermal noise both in the opened-loop 
driving and in the closed-loop driving, we find that the closed-loop driving does not reduce the RMS 
noise amplitude. We observe that the RMS noise frequency can be reduced by increasing the quality 
factor and drive amplitude in the closed-loop driving system. The experiment and simulation validate 
the feasibility of closed-loop driving and confirm the validity of the averaged equation. The 
experiment and simulation results indicate the electrical noise of closed-loop driving circuitry is bigger 
than the mechanical-thermal noise and with the driving mass decreasing, the mechanical-thermal noise 
may get bigger than the electrical noise of closed-loop driving circuitry. 
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